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We summarize the results for the master integrals of the three-loop quark and gluon form factor 
in massless QCD. Working in dimensional regularization we extract poles up to 1/ e 6 . The com- 
putational techniques involve, among others, the expansion of higher transcendental functions 
and the Mellin-Barnes method. The coefficients of the Laurent expansion in e are given either 
analytically or numerically to high precision. 
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1. Introduction and computational methods 

The quark form factor y* — s> qq and gluon form factor H — s> gg (effective coupling) are the 
simplest processes containing IR divergences at higher orders in massless QFT, and therefore are of 
particular interest in many aspects. They have been used to predict the IR pole structure of multi-leg 
amplitudes [1-4]. The form factors can also be exploited to extract resummation coefficients [5,6], 
and they enter the purely virtual corrections to a number of collider reactions (Drell-Yan process, 
Higgs production and decay, DIS). Besides phenomenological applications, a major motivation for 
obtaining analytic results at three-loop order and beyond is finding and understanding structures in 
massless gauge theories that generalize to an arbitrary number of loops. Much progress has been 
achieved in the prediction of all-order singularity structures in QCD [7-9], in conjectures about the 
all-orders behaviour of maximally supersymmetric Yang-Mills theories (see e.g. [10-12]) and in 
investigations of the finiteness of N=8 supergravity, see e.g. [13-15]. 

The two-loop corrections to the massless-quark [16-18] and gluon [19,20] form factors were 
computed in dimensional regularisation with D = 4 — 2s to order e° and subsequently extended 
to all orders in e in ref. [21]. The three-loop form factors to order £ _1 (and £° for contributions 
involving fermion loops in the quark form factor) were computed in [6,22]. Recently, also the 
three-loop form factors through £° became available [23] (see also [24]). 

In order to calculate the quark and gluon form factors at higher orders in perturbation theory, 
the amplitudes are reduced to a small set of master integrals by means of algebraic reduction pro- 
cedures [25-31]. At the three-loop level, the reduction results in 22 master integrals. Eight of them 
are products of one-loop and two-loop vertex functions or three-loop two-point functions, both of 
which are known to sufficiently high orders in e [21,25,26,32,33]. The remaining fourteen masters 
are genuine three-loop vertex functions which are depicted in Fig. |l|. Each topology contains only 
one master integral, and corresponds to two-particle cuts of the master integrals for massless four- 
loop off-shell propagator integrals [34]. Working in dimensional regularisation with D = 4 — 2e 
and expanding the master integrals in a Laurent series in e, the finite part of the three-loop form fac- 
tors requires the extraction of all coefficients through (polylogarithmic) weight six, i.e. coefficients 
containing terms up to n 6 or Q. 

The computational methods that we use during the calculation are manifold. The easier inte- 
grals (A5 1, A5 2, Ag.i) contain only T-functions and their expansion is straightforward. The more 
complicated masters Ag3, A 7 j, A 7 2 , and A 7 4 are also represented in a closed form in terms of 
hypergeometric functions of unit argument. The latter are expanded in e by means of HypExp [35, 
36]. The remaining integrals possess multiple Mellin-Barnes (MB) [37-41] representations. Their 
analytic continuation to e = was done with MB [42] and MBresolve [43]. These packages were 
also used for numerical cross checks. In addition we performed numerical checks with the sector 
decomposition methods of [44,45] and the FIESTA [46] package. 

As analytic techniques we apply Barnes's lemmas and the theorem of residues to the multiple 
MB integrals, and insert integral representations of higher transcendental functions where appro- 
priate. We also make use of the HPL [47] and barnesroutines [48] packages, as well as the 
nested sums algorithm [49,50]. 
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2. Results 

Below, we list the results [51-53] for the fourteen genuine vertex-type master integrals from 
Fig. g} We use the following definitions, 

and moreover we tacitly assume that all propagators contain an infinitesimal +ir] with r\ > 0. 

2.1 Integrals with five or six propagators 

We start with three integrals that can be displayed in a closed form in terms of T-functions only. 
Their expansion about e = can be easily performed with standard computer algebra programs. 
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The next integral, A6,2» can be written in terms of a two-fold Mellin-Barnes representation. 

fdwi fdw 

r(l-2e)r(2-4e) JlniJlni 
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The contour integrals in the complex plane can be chosen as straight lines parallel to the imaginary 
axis, i.e. the real parts c\ and c 2 along the curves are constant. Their values, together with that of 
6, must be chosen such as to separate left poles of T-functions from right ones [38-40], which is 
achieved by taking c\ = —6/5, c 2 = —1/2, and —1/15 < 6 < 3/20. The analytic continuation to 
6 = is therefore trivial, and after summing the residues one obtains 
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In Ref. [36], two more orders of the £-expansion can be found. The next integral, A^, can be 
displayed in a closed form valid to all orders in £ in terms of hypergeometric functions of unit 
argument. The latter are expanded in e by means of the Mathematica package HypExp [35,36], 
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2.2 Integrals with seven or eight propagators 

Since the explicit expressions for the integrals that follow are very lengthy, we summarize 
in this paragraph the respective techniques and subsequently only give the results of the Laurent 
expansion in £. Integrals Aj^, A7.2, and A7.4 reveal a closed form in terms of hypergeometric 
functions of unit argument, whose £ expansions are carried out with HypExp [35, 36]. Integrals 
A 7 ,3, A7 5, and A§ possess multiple Mellin-Barnes representations which are three-, four- and four- 
fold respectively. Only in the case of Ag the analytic continuation to £ = is non-trivial and is 
carried out with MB [42]. The explicit results read 
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2.3 Integrals with nine propagators 
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Integrals A9 1, A9 2, and A9 4 possess a six-fold MB representation each. The analytic con- 
tinuation, carried out with MB [42], involves in each case approximately 200 steps. We apply the 
techniques described in section [l], and stress that all results were obtained by purely analytic steps. 
We have 
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The numbers were obtained with MB . m [42]. We stress that all other terms in fl2.14| ) - ( |2. 16| ) were 
derived by purely analytic steps. The analytic result of the simple pole of A9 4 can be extracted 
from [23]. It turns out that for each of the above integrals a corresponding one with an irreducible 
scalar product in the numerator can be chosen whose coefficients of the £ -expansion have homo- 
geneous weight. The suitably chosen numerators are r 2 for A91 and A9 4, and (I — p\) 2 for A9 2. 
The property of homogeneous weight is very helpful when one uses the PSLQ algorithm [54], by 
means of which we obtain 
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Since each topology contains only one master integral, each A9 j is related to its corresponding 

In) 

A 9 {. We established these relations with the Laporta algorithm [27, 28, 30] and checked that they 
are filfilled by the above expressions. 

The analytic expressions for the remaining coefficients are within reach. In the finite part of 
Ag,2 we are missing only ~ 30 terms, all four- and five-fold MB integrals. From what we can 
judge it is not possible to process these MB kernels by purely analytic steps. We therefore plan 
to break them down to lower-dimensional integrals over ordinary Feynman parameters and then 
use the PSLQ algorithm. In the case of A9 4 we are left with <^(10 3 ) MB terms which are at most 
three-fold. We are confident that they can be processed by purely analytic steps. 
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